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ABSTRACT 

Context. The rotational lines of carbon monoxide and the fine structure lines of ClI and of the most abundant metals, emitted during 
the epoch of enhanced star formation in the universe, are redshifted in the frequency channels where the present-day and future CMB 
experiments are sensitive. 

Aims. We estimate the contribution to the CMB angular power spectrum by the emission in such lines in merging star-forming 
galaxies. 

Methods. We used the Lacey-Cole approach to characterize the distribution of the merging halos, together with a parametrization for 
the star formation rate in each of them. Using observational data from a sample of local, low-redshift, and high-redshift objects, we 
calibrated the luminosity in each line as a function of the star formation rate. 

Results. We show that the correlation term arising from CO line emission is a significant source of foreground for CMB in a broad 
range of frequencies (in particular in the 20 - 60 GHz band) and for 1000 < / < 8000, corresponding to angular scales smaller than 
10 arcminutes. Moreover, we demonstrate that observing with different spectral resolutions will give the possibility of increasing 
the amplitude of the signal up to two orders of magnitude in C/ and will help separate the line contribution from practically all 
other foreground sources and from the primary fluctuations themselves, since these show no significant dependence on the spectral 
resolution. 

Conclusions. We propose to perform observations with varying spectral bandwidths (10"^^ < Av/vobs < 10"') as a new tool to 
construct a tomography of the universe, by probing different redshift slices with varying thickness. This should yield new constraints 
on the regions responsible for the metal enrichment in the universe and on their clustering pattern and will lead to new hints about the 
reionization epoch and the cosmological parameters, including erg. 

Key words. Cosmology: cosmic microwave background - Cosmology: theory - Galaxies: intergalactic medium 



00 

o 

O 
(N 



<N 

S: 

9^ 
6 

(N 
> 

(N 

in 

o ' 

1. Introduction from other foregrounds generated by continuum emission (dust, 

radio sources, SZ effects in clusters of galaxies, and the pri- 
The nucleosynthesis in the first generation of stars is mordial fluctuations themselves) does not depend on the spec- 
■ able to produce significant amounts of carbon and oxygen tj-^i resolution of the observations. Observations of the power 
k> ; (|Heger&Wooslev||2002t [Yoshida et al.||2006|). This bring our spectrum of CO line emission with high spectral resolution al- 
, attention to the carbon monoxide rotational lines originating narrow slices of the universe to be probed at different red- 
' from the star-forming halos in the early universe. Recent ob- shifts, with a thickness Az/z ~ Av/vobs- Performing our compu- 
servations show that these hnes are very bright in the most tations for different specti-al resolutions and proving the strength 
distant quasars and radiogalaxies at redshift z ~ 4 (see e.g. jhis method, we recognized that an analogous approach is 
iGreveetal. 2005') and also in the star-forming galaxies in our ^^ed in demonstrating the results of numerical simulations of 
vicinity (Gao et al.!|2001t IWeiB et al.| 2005b; Bayet etalj|2006t the large-scale structure of the universe. In that case the most ef- 
iBaaneLal. 2008). They_are also the brightest radio lines in our festive contrast is achieved using slices with a very small thick- 
Galaxy iWright et al. 199lt| Fixsen et al, .1999j. ness (~ 15 Mpc) compai-ed to the huge dimensions of the box 
In this paper we are making an attempt to compute the angu- (~ GpcQ)- The observations in the CO lines can use comparable 
lar power spectrum of the foreground fluctuations due to emis- Av/vobs to increase the amplitude of the signal and to separate it 
sion in such lines from merging star-forming galaxies. The main from other continuum sources. 

difficulty of this task is related to the presence of the contin- n is rather difficult, at present, to obtain a reliable theo- 
uum emission of dust at the same frequencies. To separate the retical estimate of the luminosity of diff'erent CO lines in the 
line contribution, we propose observing in several spectral bands star-forming halos (Silk&Spaans 1997; Combes et al. 199^ 
with resolutions in the range Av/vobs = 10 ' - 10 ^ We demon- iGreve & Sommer-Larsenl l2008), and for this reason we use ex- 
strate that with this technique it will be possible to increase jsting observational data to calibrate our model. We consider lo- 
the amplitude of the signal due to the lines emission signif- cal merging galaxies and luminous infrared galaxies, together 
icandy (over one order of magnitude), while the contribution with a sample of high-redshift, sub-millimeter galaxies, ob- 
served at z ~ 2 - 4. Using these samples, we calibrate the lumi- 
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nosity of the different lines on the star formation rate of the ob- 
ject, assuming a linear scaling. We follow the same approach as 
outlined in our previous paper (Righi et al. 2008), where we pre- 
sented a model for computing the fluctuations due to dust emis- 
sion in merging star-forming halos. This model, based on th e 
extended Press-Schechter formalism (EPS. lLacev & Colel 19931) . 
allowed us to obtain a statistical description of the distribution of 
the halos as a function of their star formation rate. With these two 
independent ingredients, we can compute the power spectrum of 
angular fluctuations arising from the emission in the CO lines, 
at any frequency and for any spectral resolution. This approach 
will allow the integral properties of the population of weak star- 
forming merging galaxies to be measured in a broad redshift 
range: by calibrating with observational data obtained from dif- 
ferent galaxy samples, we can describe the history of formation 
of giant molecular clouds and model the subsequent process of 
CO enrichment of the IGM throughout cosmic history. 

The beauty and the strength of the proposed method of 
observing the angular fluctuations in the narrow hne emission 
is connected to the fact that, for a given line at a fixed ob- 
serving frequency, we detect the radiation from a slice of the 
universe with given Az. The data on the cosmic star formation 
history (the so-called M adau plot, iMadau et al] 1 19961 Il998t 
iHopkins & Beacomll2006l) shows that the bulk of the star forma- 
tion activity in the universe takes place in a very broad redshift 
range. Here we are considering the star-forming objects in the 
range ^ z ^ 10, corresponding to 10 - 115 GHz for the 
first CO transitio n. This is close to the WMAP and PLANCK 
LFFs detectors jBennett et all 120031: iBersanelli & Mandolesl 
12000 *). Future ground-based experiments operating in the 
centimeter band and able to perform observations in many 
relatively narrow spectral channels (10^^ < Av/vobs < 10"') 
might open the possibility of separating the contribution of line 
emission from any other continuum source using the method 
propose d in this p aper Mo re accurate s emi-analytical models 

(see e.g.lGranato et al. 2004tlLacev et al.ll2008: Cole et al 2008.: 

iReed et al.ll20'08ir might allow the predictions to be extended to 
even higher redshifts and to correspondingly lower frequencies. 

Higher CO transitions are still bright up to ~ 350 GHz 
and are therefore in teresting for the PLANCK H FFs bands 
(Lamarr e et al. 2003h. as well as for ACT and SPT (iKosowskvl 
r2006tiRuhl et al.ll2004 . The same spectral band can cover dif- 
ferent slices of the universe, and from this point of view we will 
see that the 10-40 GHz region is the most interesting one for the 
CO emission, because only the first one or two transitions (and 
the corresponding redshift slices) contribute significantly to the 
C/'s in this band. We present, however, the frequency (and red- 
shift) distribution for all the CO transitions up to 7 - 7, which 
will permit recognition of the redshift range that gives the main 
contribution to the signal. 

According to the same model, we also compute the dififeren- 
tial source counts curve for the emission in every line. The cumu- 
lated signal from many weak sources might contribute more to 
the power spectrum due to the clustering. Therefore, in addition 
to the simple Poisson fluctuations, we computed the correlation 
signal taking into account the regions of higher merging activ- 
ity, corresponding to the positions where the future clusters and 
super-clusters of galaxies will appear. 

From the same sample of objects that we used to calibrate 
our model, we see that the emission from the CN, HCN, HNC, 
and HCO^ molecules, with resonant frequency close to the CO 
ones (within ~ 20%), d o not add more t han 10% to the am- 
plitude of the CO signal (iBaan et al.ll2008l) . At higher frequen- 



cies, the dominant contribution comes from the CII 158 fim line, 
while other atomic lines are much weaker Here, however, we 
are interested not in the absolute luminosity of the lines, but 
in their ratio to the primordial fluctuation and to dust contin- 
uum emission for the same angular scales, in the same spectral 
band and with a given Av/ Vobs • This argument makes the first two 
CO transitions more important than the (more luminous) higher 
transitions, given their contribution to the power spectrum of an- 
gular fluctuations relative to the primordial CMB signal and to 
the foreground from extragalactic sources emitting in the sub- 
millimeter band due to the presence of large amounts of dust. 
There are no galactic foregrounds that can mimic a similar de- 
pendence on the spectral resolution to the signal from CO. The 
first CO transitions are usually saturated, and their ratio follows 
the same Rayleigh- Jeans law of the CMB at low frequencies. At 
the same time, the dust particles have very low opacity at low fre- 
quencies making their effect weaker with decreasing frequency, 
compared with CO (1-0) and (2-1) transitions. 

We use the EPS formalism because it is relatively simple 
and permits us to demonstrate the key results of the proposed 
method. The deviations that the EPS intr oduces with respe ct to 
the predictions of numerical simulation jCole et al.l l2008l) are 
still smaller than the uncertainties in the observational sample 
that we use for the calibration of the CO line luminosities. Such 
uncertainties are mainly due to the data for the star formation 
rates: in the low-redshift sample we use only the IRAS data, 
combined with the Kennicutt relation, to estimate the star forma- 
tion rate. The galaxies in the high-redshift sample, on the other 
hand, often harbor an AGN and this is a source of further uncer- 
tainty. An additional difficulty comes from the different angular 
resolution of the CO and infrared observations. 

Future observations with different bandwidths will measure 
the dependence of the CO contribution on spectral resolution. 
More precise numerical models can refine the predictions of our 
simple model. In our view, it is important to propose a new 
method of directly measuring the rate of enrichment of the uni- 
verse with CO not in the very rare and extremely bright objects 
(as QSO and radiogalaxies at z ~ 4 - 6), but in the most abundant 
and less bright objects, which cannot be detected individually. 



2. Star formation model 

Following the same strategy as in our previous work on the in- 
frared sources (Righi et al. 2008), we model the star formation in 
a merger episode and charac terize the d istribution of the merging 
using the EPS formalism CLacev & Col&,1993.) . In a very recent 
paper. iNeistein & Dekell (l2008l) improve the EPS formalism, but 
their new estimates of the merger rates only deviates signifi- 
cantly from the results of Lacey & Cole for the minor mergers. 
For the merging mass ratios considered in our approach (0. 1-10), 
the discrepancy in the two models is only 20%. In the following 
we briefly outline the main equation of our star formation model. 

In each merging episode, a given amount of gas is converted 
into stars. Considering two halos of mass Mi and M2, merging 
to yield a halo of mass M = Mi -1- M2, a given amount of the halo 
mass is converted into stellar mass 



I Qb M2 

Mi 77 Mi . 

* Qm ' M/2 



(1) 



The amount of stellar mass formed from the halo Mi is therefore 
proportional to its mass and to the mass of the merging halo M2, 
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Fig. 1. The distribution of the % - LjM ratios as a function of the upper rotational quantum number of the CO transitions (the 
corresponding frequency is shown in the upper axis in GHz). Three different samples are considered: low-re dshift IRAS galax- 
ies (blue points), high-redshift sub-millimeter galaxies (red points), and high-redshift radiogalaxies and QSOs ("Greve et al.''2005|; 
ISolomon & Vanden Boutl2005l green points). M82 and the Antennae are shown in black. Our Galaxy ( Wright et al. 1991; Cox 200^ 
is represented with the cyan triangles. See Tables [1] and |2] for details. The black solid line is SED for M82, for which we have data 
for the full set of lines. Dashed lines represent the upper and lower limits for % in the sample. 



if this is sufficiently massive. In the same way, for the second The total stellar mass produced in this merging episode is then 
halo, the sum 

M.^Ml.Ml.4^j^. (3) 

The parameter rj is the star formation efficiency. We take this 
2 fib Ml to be 5%, in order to match the observations of the cosmic star 

^* = '? ^2 ■ (2) formation history ( 'Hopkins & Beacomll2006l) . 
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Table 1. Luminosities and H ratios (in boldface) of the CO lines for the local and low-redshift sample of galaxies. 



Object name 


Type 




SFR 


CO (1-0) 


CO (2-1) 


CO (3-2) 


CO (4-3) 


CO (5-4) 


CO (6-5) 


CO (7-6) 










115.3 GHz 


230.5 GHz 


345.8 GHz 


461.0 GHz 


576.3 GHz 


691.5 GHz 


806.7 GHz 






[ivipcj 




\j ^ 
l^oJ 


\T 1 

I^-gJ 












M82 


SB 


3.63 


10.1 


3.7 X 10^ a 


2.8 X 10= a 


7.1 X HP a 


9.8 X UP a 


9.7 X 10= a 


9.6 X 10= a 


9.0 X 10= a 










J. / X lU 


Z.o X lU 


/.U X lU 


O 7 v Ift4 
v. / X lU 


a V ift'i 
y.o X lu 


Q C V 1 ft"* 


fi O s/ 1 ft4 






M82 Ldusi 


J Leo ratio: 




0.005 


0.009 


0.02 


0.05 


0.12 


0.24 


Antennae 


SB 


21 8 


118 


1 S V UP K 
1 .o X 1 U D 

1.5 X 10* 














Milky ^Vay 
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O.J X lu C 

2.8 X lO"* 


1 2 X 10^ c 
4.0 X 10" 


1 3 X 10^ c 
4.3 X 10" 


9 5 X 10"^ c 
3.2 X 10" 






Ti? A c n 1 (Yin 1 ncvi 
IK A J UiU / / -I /u/ 




1 /I «; 


/Z.V 


z.y X lu a 
4.0 X 10' 














IKAj U1j04-1U4Z 




zu / 


QC /I 


Z.J X lu a 
2.3 X 10' 














TDACn/i/i«;/i /icjc 




Ton 


IUj 


1 .0 X iu a 
1.5 X 10' 














TD A c nc'con AC'cn 

IKAj UojzU-OojU 




ZUj 




I . J X 1 u a 
1.3 X 10' 














iKAo UV 111-1 UU / 




243 


171 


1 s V in5 

J. a X lU a 

2.2 X 10' 
















CK 


ICS 


j41 


AOs/ 1 n5 A 
o.v X lu a 

Z.U X lU 














IRAS 14378-3651 




315 


242 


2.7 X 10' d 
1.1 X 10' 














IRAS 18293-3413 




80.6 


110 


5.3 X 10' d 
4.8 X 10' 


2.2 X 10" d 
2.0 X 10" 












IRAS 19115-2124 




216 


127 


5.1 X 10' d 
4.0 X 10' 














IRAS 20550+1656 




136 


127 


1.7 X 10' d 
1.3 X 10' 


4.8 X 10' d 
3.8 X 10' 












IRAS 22491-1808 


SB 


350 


220 


3.0 X 10' d 















1.4 X 10' 



References: (a^l Weifi et all I2005hl). fMIGao et all 1200 Ih. (c)IWright et anil991h. (d)lBaan et all 120081). 



The lifetime of the starburst phase is then introduced 
following the results of th e several numerical simulations 
available in Uterature (e.g . iMihos & Hernauistj [l994l Il996t 
ISpringel & HernquistI l2005l iRobertson et all l2006l) . There are 
two star formation bursts, corresponding to the first close pas- 
sage of one galaxy around the other and to the final coalescence, 
with a typical duration of ~ 300 and ^ 100 Myr, respectively. 
A third phase, identified as the star formation activity already 
present in the two galaxy before the interaction, is not included 
in the model, since it is not directly trigger by the merging. This 
simple parametrization allows us to derive a statistical descrip- 
tion of the merging events as a function of their star formation 
rates. 

Comparing the typical timescale of the bright phase (~ 100 
Myr) with the characteristic evolution time for old stars, which is 
about 1 Myr, one can easily prove that these objects are rapidly 
enriched with metals already at high redshift. This is the main 
argument that motivates the work presented in this paper. 

The next step is then to connect the star formation rate with 
the luminosity in a given band. In our previous work, we were in- 
terested in describing a population of far-infrared sources, so we 
used the Kennicutt relation for that band (Kennicutt 1998). Here 
we aim to describe the luminosity in the line for a set of molecu- 
lar and atomic transitions that occur in the observing band of the 
CMB experiments. We assume that line luminosity scales with 
the star formation rate 

Aine[io] =^-M*[Mo/yr]. (4) 

The constant 'R will be calibrated in the next section, taking a 
sample of objects of different nature into account both at low and 
high redshifts, along with their observed value of the luminosity 
in each line. 



3. Line luminosities from a sample of objects 

In the previous section we introduced a simple relation between 
the star formation rate and the luminosity in a given line, sim- 
ilar to the existing relations in different bands (see iKennicuttI 
1998, for a review). The constant factor connecting the two quan- 
tities depends on several parameters, and it is not easy to deter- 
mine it. To get a rough estimate of its value for different lines, 
we consider a sample of objects and retrieve data on their line 
luminosity from the literature. 

For the CO transitions, we consider both a low- a nd a high- 
redshift sample. The first is derived from the list of iBaan et al.l 
(2008), which contains more than a hundred sources, observed 
with different instruments. We selected in this list only the galax- 
ies with optical size smaller than the beam size of the observa- 
tion. This is to avoid an underestimate of the flux due to incom- 
plete coverage of the source. More detailed studies have been 
done on well-known local objects like M82 (IWeiB et al. 2005b|) 
and NGC4038/39 (also known as The Antennae, Gao et all 
l200Tt iBavet et^^M^ . These are, however, extended sources 
so one has to be very careful with the observational results, since 
they depend on t he actual region co vered by the instrument. In 
the case of M82. IWeiB et al.l(l2005bh present an accurate model 
of the CO emission from this galaxy, which takes both the cen- 
tral and the outer regions into account. The combined emission 
is a reliable estimate of what would be observed if M82 was 
shifted to cosmological distances, therefore we will use it to 
calibrate the ^ ratio in our model. Th e relatively low mass of 
M82 (M ~ 10'° Mn. i Mavva eTai1l2006 ) is typical of the objects 
that contribute to the bulk of the correlation term estimated here 
(M = 10'^ - 10" Mo for redshift z = 3 - 5). The data of the 
low-redshift sample are presen ted in Table jT] the star formation 
rates are obtained applying the iKennicufl ( Il998l) relation to the 
far-infrared (8 - lOOOyum) luminosity given in the IRAS revised 
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Object name 


Type 




SFR 


CO (1-0) 


CO (2-1) 


CO (3-2) 


CO (4-3) 


CO (5-4) 


CO (6-5) 


CO (7-6) 










1 i J.J urlZ 


ZjU.J OrlZ 


J4J.O OrlZ 


zlAl A ni-Iv 
40 i .\i urlZ 


J /O.J yjiiL 


Oy 1 .J OrlZ 


SHA 1 ni-Iv 








[Mo/yr] 




[i-ol 


[i-ol 








[i-o] 






Z. jZ 


500 a 






o.u X lu D-a 
1.2 X lO" 


1.1 X 10' b 

2.2 X 10" 


1 .7 X 10' b 
3.4 X 10" 


1 A s/ 1 A' K 
1 .0 X lU D 

3.2 X 10" 


1 /I V 1 A' K 
1 .4 X lU D 

2.8 X 10" 


ciMM JUzJyo* 


AGN 


1.06 


975 e 




"7 A ^ 1 a6 C 

/.9 X 10 r 
8.1 X 10^ 












CIVfllVfl 7 111 OA 


AGN 


3.41 


810 g 


1 A 1 a7 ~ 
l.U X lU g 

1.2 X lO-* 






1 "7 ■^y' 1 A8 f 

1 . / X lu r 
2.1 X 10' 








T 1 (^"Xf^f^ 
jIVlIVl J lOJOO 




z.4j 


14jj 1 






"7 7 1 a7 f t, 

/. / X lu r,n 
5.3 X 10" 








1 1 1 A8 K 

J.J X lU n 

2.3 X 10' 


jiVliVl J lOJ / 1 




Z.JO 


877 i 






4.U X lu r 
4.6 X 10" 










ciMM Jzzl /4 


SB 


3.10 


1800 e 






5.1 X lO' f 
2.8 X 10" 










jMM JU44J 1 * 


C D 1 A M 

CiD + AOIN 


2.51 


450 e 






1 /I \y 1 a7 t, ; 

1 .4 X lu n,j 
3.1 X 10" 








/I V 1 a7 V, 

Z.4 X lu n 
5.3 X 10" 


ciMM JUy4il* 


SB + AON 


3.35 


1200 g 


< Z.J X 10 g 

< 2.1 X 10^ 






o A ,y 1 a7 t. ; 
0.9 X 10 n,J 

7.4 X 10-* 








jMM J lOJOo 


Sd + AOIN 


2.38 


897 i 






A "3 \y 1 a7 t, ; 

y.j X lU n,j 
1.0 X 10 -^ 








J. / X lU n,j 
4.1 X 10' 




CD _i_ A (""M 


Z.oU 


'inn v 






A A v 1 A^ l- 1 
o.o X lU K,l 

1.3 X 10' 










T1 A(\\ 1 J, 




Z.JO 


jdU e 






Z.4 X lU 1,111 

6.7 X 10" 








A O S/ 1 A^ m 

o.y X lu m 
1.9 X 10' 


T1 O'^^AQ 


JO + A01>l 


2 20 


1 163 n 






J .o X 1 u n 
4.8 X 10" 






1 A V 1 n8 Vi 

1 .o X iu n 
1.4 X 10' 




CDi^l T 1 (^A 




1 .44 


1539 o 


J.Z X lU 

2.1 X 10^ 


1 . J X lu a 
9.7 X 10^ 






J.o X lU p 

2.3 X 10" 






LiUUUj J1zjoj4 


SB 


1 .22 


950 q 




Z.D X lU q 

2.8 X 10* 












IjIZ* 


LBG 


2.73 


15 e 






j.y X lu r 
4.0 X 10" 










c\ no*; 7 J. 
(J UV J / * 




1 /1 1 

1 .41 


vuu e 




"7 /; w 1 a6 _ 
/.O X lU s 

8.5 X 10' 














QSO 


2.29 


540 e 






Z.Z X lU t-W 

4.0 X 10" 






4.Z X lU V 

7,7x10" 








Z. JD 


o lu e 






4.0 X lU x-da 

5.9 X 10" 


1 V 1 A^ 1 
1 .Z X lU Z 

1.5 X 10' 


1 "7 w 1 A^ -F 
1 . / X 1 1.) Z 

2.1 X 10' 




4.0 X lU Z 

5.7 X 10' 


D V TAA 1 1 _i_ 

K A J Uv 1 1 ■* 


QSO 


2.80 


345 e 






/.I X lu ao 
2.1 X 10" 














Z.o4 


jouu e 






Z.J X lu ao 
6.4 X 10" 










AT^ 1 _1_ 


QSO 


3.20 


435 e 








J.Z X lu ac 
7.2 X 10" 








PSS J2322* 


QSO 


4.11 


1800 e 


2.6 X 10"^ ad 
1.4 X 10' 


2.5 X 10' ad 
1.4 X 10" 


2.3 X 10" ae 
1.3 X 10' 




2.5 X 10* ae 
1.4 X 10' 






BRI 0952* 


QSO 


4.43 


360 e 










4.3 X 10' af 
1.5 X 10' 






B3 J2330 


HzRG 


3.09 


1950 e 








1.1 X 10* ag 
5.6 X 10" 








TN JO 121 


HzRG 


3.52 


1050 e 








1.3 X 10* ah 
1.2 X 10' 








6C 1909 
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bright galaxies sample (ISanders et alJl2003l) . Distances are from 
the same catalog: conversion from proper to luminosity distance 
has been computed using the cosmology-corrected redshift in 
the Nasa Extragalactic Database (NED). For our Galaxy we use 
the value of star formation rate given bv [Coxl(l2000l) . 

In the same table, we also present the ratio between dust 
and line luminositie s in M82, using the spectru r n given in 
iLagache et all (|2005b and the mo del of l^jBeLaLl (l2005bh . In 
our previous paper (Righi et al.ll2008h . we demonstrated that 
dust emission from merging galaxies is the most important high- 
frequency extragalactic foreground, therefore it is very important 
to compare with it. As we can see from the corresponding line 



in Table[Tl the ratio of dust-to-CO-line emission grows with fre- 
quency; therefore, the first CO transitions are more effective and 
should be favored from the point of view of their ratio to the dust 
continuum. 

The high-redshift sample contains many sub-millimeter 
galaxies. These are very massive and luminous systems, pow- 
ered by intense starburst and/or AGN, and they emit the bulk 
of their radiation in the infrared band. This activity is triggered 
by merging episodes that provide large amount of gas and re- 
sults in high values of the star formation rate. We retrieved from 
the literature the line luminosity for several objects, including 
QSOs and radiogalaxies. These are reported in Table |2l lumi- 
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nosities of the lensed sources (marke d with ★) are corre cted 
with the magnification factors given in iGreve et aTl (1200 5*) and 
[Solomon & Vanden Bout (2005). Where more than one obser- 
vation were available we assumed an average between them. 
These observations were carried out with interferometric tech- 
nique; therefore, they are more sensitive to bright small spots 
and not to the continuum emission. Comparing the average 
of this sample with the one obtained for M82, we see that there 
is no big difference, though the CO spectral energy distribution 
for the high-redshift objects tends to peak at higher values of the 
rotational quantum number J. This may derive from differences 
in the temperature and in the density of the gas. Furthermore, in 
the high redshift objects, the presence of the non-negligible ra- 
diation field of the CMB contributes to increasing the excitation 
level of the CO molecule. 

We summarize all these observational data in Fig. [1] where 
we plot the H - L/M ratio as a function of the upper rotational 
quantum number of the transition. A higher value of 'R corre- 
sponds to a higher luminosity of the galaxy per unit star for- 
mation rate, hence to a higher amplitude in the fluctuations. We 
present the two samples of objects described above: local ob- 
jects, low-redshift IRAS galaxies, high-redshift sub-millimeter 
galaxies, QSO and radiogalaxies. The spectral energy distribu- 
tion of M82, for which we have the full set of lines, is overplot- 
ted. We show here also the values for our Galaxy, as observed 
by COBE/FIRAS (Wright et al. 1991), assuming a star forma- 
tion rate of 3 M^/yr (ICoxll2000l) . From this dataset we extract 
an estimate of the lower and upper value of fi in every transi- 
tion. Notice that local merging objects, which are intrinsically 
weaker than the QSO, are close to the upper limit, meaning that 
they have a higher luminosity per unit star formation rate. M82 
and the Antennae, for example, are extremely bright in the first 
two transitions, compared with the other objects. There is rather 
large scatter in the distribution of these points at different J. This 
might be due to the uncertainties in the CO properties or more 
probably in the value of M measured in these objects. 

The observations of bright distant quasars and galaxies show 
that they often have carbon and oxygen abundance close to or 
even higher than the solar value (Fan 2006). It is clear that to 
the average abundance of metals is much lower in the high- 
redshift universe than in the regions with intense star formation, 
where it can increase extremely rapidly. The Magellanic clouds 
(LMC and SMC) have significantly lower chemical abundance 
than our Galaxy; nevertheless, the luminosities of the CO line s 
are comparable dCohen et al.1 119881; ICaldweU & Kutnedll996h . 
This problem is of particular interest for the first line of the CO, 
since it is heavily saturated and therefore its brightness depends 
weakly on chemical abundance. On the other hand, higher-7 
lines are usually optically thin and therefore their brightness 
strongly depends on the chemical abundance. This dependence 
might be non-linear, with an index close to 1 .5 - 2, leading to 
lower CO emission in metal-poor systems. This is also con- 
nected with the low dust abundance and a consequently less ef- 
ficient shielding from the ultraviolet radiation, which destroys 
CO molecules (Combes, priv. comm.). However, in this paper 
we are dealing with star forming regions, in which rapidly evolv- 
ing massive stars efficiently enrich the medium with significant 
amounts of metals. M82, which is the object used to calibrate 
most of the results presente d here, has a metallicity 1 .5 times the 
solar value dGalUano et al.l l2008). 

The saturation of the first two CO lines might lead to a non-linear 
dependence of the li ne lumi nosity wit h the star formation rate. 
This was shown by iNaravanan et al.l (|2007|) : using a radiative 
transfer model combined with numerical simulations they found 
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Fig. 2. The actual comoving width of the slice (in Mpc) as a 
function of redshift, probed by different values of the spectral 
resolution Av/vobs- 

that this dependence is linear for the higher transitions {J > 3), 
but might become non linear at lower J. 

4. The angular power spectrum 

In the Appendix |A] we present a detailed derivation of the ex- 
pressions we used for the angular power spectrum generated by 
the emission in these lines. Here we briefly outline the main re- 
sults. First of all, we should note that the measured signal /v(n) 
is actually a convolution of the underlying signal Iy(r = rn) with 
an experimental window function S(r): 

/v(r) = J flfr'S(r-r')/v(r')- (5) 

This window function S simply accounts for the finite angular 
and spectral resolutions of the observing device. The variance of 
the signal can be expressed as a function of the measured power 
spectrum (|/^kft 

^'"^^/ (^l^^.kl'=/(^l^kl^^..«l^^l'' 

where ly denotes the volume weighted intensity amplitude for 
a given source, and Ps(k) is the power spectrum describing the 
spatial distribution of the sources. This integral reflects the con- 
tribution of anisotropy (given by the power spectrum Ps(k)) to 
each k scale range, and this depends on the size and distribu- 
tion of the sources. If these are randomly (Poisson) distributed in 
space, then Ps{k) - 1/n, with n the average source number den- 
sity. This means that all scales contribute with the same amount 
of anisotropy (as given by the A:-independent power spectrum 
Pj), down to a minimal linear scale (maximum k) correspond- 
ing to either the size of the source or the resolution element of 
the instrument. If the source size is taken as arbitrarily small, 
and we fix the angular resolution, then the variance or measured 
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Fig. 3. The predicted differential source counts for the CO lines at 30, 70, and 100 GHz, for a spectral resolution Av/vobs = 10"^. The 
vertical line is the expected sensitivity of ALMA to the line emission, computed for the same spectral resolution and an integration 
time of 3 hours. 
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Fig. 4. The contribution to the cosmic microwave background 
radiation (in temperature units) from the different CO lines and 
of their sum (black line), compared with the signal from dust 
(orange). 



scale, then further improving such resolution will make no dif- 
ference. This can be rephrased in Fourier space as follows. If 
sources are distributed such that their power spectrum is propor- 
tional to that of the underlying density field (Psik) oc P„(k)), 
then for some large k the contribution of the power spectrum 
to the integral of Eq. ^ will be negligible (since for large k 
we have that P,„(k) — > 0). In this scenario, increasing the spec- 
tral resolution will not change the integral over Ps(k). These two 
distinct regimes will be addressed in detail below. 

After projecting the power spectrum on the sphere, we find 
expressions for the Poisson (Cf) and correlation (Cf) angular 
power spectra (see Appendix). We recall that their relation to the 
angular correlation function reads as 



<A,(ni)/v(n2)) 



J]^{Ci^Cf)PKn,-n2), 



(7) 



with P/(ni ■ n2) the Legendre polynomia of order /. In this ap- 
proach the resolution is expressed in terms of a comoving width 
Ar probed at a given redshift by a certain Av/vobs- It is straight- 
forward to derive the relation 
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where Av/vobs refers to the observing frequency and E{z) - 
[ilm(l + z)^ + Qa]'^^. From this equation, it is clear that, for 
a fixed spectral resolution, the actual comoving width decreases 
with increasing redshift, as it is shown in Fig.|2] 



anisotropy will increase if the instrument is sensitive to smaller 
and smaller scales in the radial direction; i.e., it will scale as 
^max ~ ^li^z) ~ l/(Av)inst, with z denoting the linear scale along 
the radial direction. More generally, if sources are Poisson dis- 
tributed, any improvement in the spectral/angular resolution that 
makes the experiment sensitive to smaller scales (that are still 
larger that the typical source size) will yield an increase in the 
measured power 

However, we know that, in the universe, sources tend to be 
clustered in regions of ~ 15-25 Mpc (comoving size) that 
eventually will become superclusters of galaxies. The distribu- 
tion of these regions will introduce more anisotropy, but only on 
scales that are actually larger than Lc- If the frequency/angular 
resolution of the experiment is at some point able to resolve this 



5. Results for the CO lines 

5.1. Source counts and background intensity 

In Fig. [3] we present the differential source counts in the CO 
lines, predicted by our model at 30, 70, and 100 GHz, for a 
spectral resolution Av/vobs - 10"^ and using M82 to calibrate 
the fi ratio. The vertical lines in each plot show the sensitivity 
of ALMA to line emiss ion (5cr value, for 3 hours integration. 
iTarenghi & WilsonI 12005.) : this value is approximately constant 
at these frequencies. The flux in the line is computed as 



'"""47rr2(l -H z) _ r , , Lho\4>th,v 

-p- - aVobsWnsti--; — y-. - 

J c/vobs^instr J 47rr^(l-Hz) 
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Fig. 5. The correlation signal for the CO emission lines at the LFI frequencies and for a spectral resolution Av/vobs = 0.2, using 
the H ratio from M82. The black line is the primordial signal of the CMB divided by 1000, the orange line is the signal from dusty 
merging star-forming galaxies. Other colors identifies the lines as indicated by the labels. 
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Fig. 6. The change in amplitude of the correlation (solid line) and 
Poisson (dashed line) signal for the CO (1-0) line for different 
values of the spectral resolution Av/vobs, at 30 GHz. The red 
and cyan lines are almost superimposed, meaning that, for these 
resolutions, the signal has already reached a constant value. The 
Poisson term, on the other hand, grows linearly with decreasing 
values of the spectral resolution. The primordial signal (black 
solid line) has been divided by 1000. 

47rr2(l +z)(Av)i„,„ 47rr2\vobs/ ' 

where Ly - L\,o\/v, Lboi is the bolometric luminosity of the line 
and r the comoving distance. This expression shows the explicit 
dependence of the measured flux on the spectral resolution of the 
instrument. The amount of sources at a given frequency is also 
dependent on the spectral resolution, just because increasing the 
value of Av/vobs corresponds to probing larger volumes of the 
universe. 



Fig. 7. The dependence of the ampUtude of the correlation signal 
on the spectral resolution of the observing instruments, for the 
CO (1-0) line at 30 GHz and for several value of the multipole 
index /. 



Given the source number counts, we can compute the contri- 
bution of each line to the cosmic microwave background temper- 
ature (Fig. nil. The total emission of the CO lines is practically 
constant with frequency in the very broad 20 - 200 GHz range. 
Its value is close to 1 fiK, and it is thousand times weaker than 
the present day uncertainties in the temperature of monopole ac- 
cording to COBE/FIRAS results (Fixsen & Mather 2002). The 
first CO line peaks at frequency ~ 15-20 GHz, corresponding 
to an emission redshift z ^ 5 - 7. The dust emission from merg- 
ing galaxies contributes two or three times more at frequencies 
around 15-20 GHz, increasing then to a value of 100 juK at 
200 GHz. 

At 30 GHz, most of the contribution to the background in the 
first CO transition is due to low-flux sources (Sy - 10""* - 10"^ 
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Fig. 8. The ratio of the correlation signal from CO emission to 
the primordial CMB signal at / = 2000 as a function of frequency 
and for a spectral resolution Av/vobs = 10"^. The black solid 
line is the sum of the correlation from the first seven CO lines, 
computed by assuming the H ratio for M82. The blue short- 
dashed line is obtained with the model based on simple Press- 
Schechter distribution. The black long-dashed line is the signal 
from CII 158 /vm line. The contribution from dust emission in 
merging galaxies is shown with the orange line. The gray point 
is the Poisson at 33 G Hz from radio sources, estimated from the 
Ide Zotti et al.l ( l2005h model, assuming a cutoff of 0.1 mJy. The 
frequency dependence of this signal (gra y hne) has been com - 
puted assuming Sv ^ v with a = 0.4 dToffolatti et al.ll2005l) . 
The shaded region around the curve for the first CO line repre- 
sents the range of uncertainties in the correlation signal, accord- 
ing to the different values of "R for the sample of objects in Fig.[T] 
and Tables [T] and |2] The first CO line is especially suitable for 
observations of higher redshifts because, at higher frequencies, 
we simultaneously observe the contributions from two or three 
slices of the universe. 



mJy), which cannot be detected directly. But since their contribu- 
tion to the correlated angular fluctuations is high, observations of 
the C/'s with different spectral resolution provide us with a pow- 
erful tool for studying the properties of the high-redshift pop- 
ulation of CO-emitting galaxies and its evolution in time. We 
discuss this in the next section. 

5.2. Power spectrum of angular fluctuations 

In this section we present the angular power spectra arising from 
the line emission of CO in the star-forming regions. Such a 
power spectrum is the sum of a Poisson term and a corre- 
lation term Cf . 

The emission lines from carbon monoxide are associated 
with the rotational transitions between different states J. Their 
emission rest frequencies lie in the range 100 - 800 GHz making 
them of particular interest for the low-frequency experiments. 
Unfortunately, there are few observations of this molecule in the 



sample that we have selected, so we were only able to retrieve 
the full data set for M82. We use this object to calibrate the H 
ratio in our computations (see Table [T]for details), but we show 
also the expected upper and lower limits in the amplitude of the 
signal, according to the value presented in Fig.[T] The correlation 
terms of the angular power spectrum for each CO line are shown 
in Fig.|5]at the three frequencies of PLANCK' s LFI instrument 
and for the corresponding spectral resolution Av/vobs - 0.2. 
As expected, the signal at low frequencies is dominated by the 
lines corresponding to the first two transitions (2-1 and 1 - 0), 
while higher states are important only at higher frequencies. This 
might appear paradoxical, since the SED of CO peaks at higher 
J. One has to bear in mind, however, that the power spectra 
are presented in temperature units and normaUzed to the CMB 
blackbody. 

Another interesting aspect is that the emission from the first 
two transitions of the CO might be slightly amplified by the con- 
ti-ibution from CN, HCN, HNC, and HCO+. These molecules 
form in the same environment of CO and has a very similar struc- 
ture. Their resonant frequencies are very close (within ~ 20%) 
to the 115 and 230 GHz transitions of CO. Their signal is then 
summed to the CO, with a slight shift in frequency and redshift, 
which leads to a smoothing effect. The luminosity o f such lines 
in nearby objects is typically ~ 10% of the CO lines dBaan et al.l 
200&) . 

As we have seen, in the literature there are also several obser- 
vations of high-redshift sub-millimeter galaxies. These objects 
are thought to be the result of very massive and gas-rich mergers, 
with very hig h star formation rates, o n the order of 10^ Mq/jt. 
According to Alexande r et alj(l20 05a b) around 40% of the sum- 
millimeter galaxies (SMG) population hosts an AGN, but these 
only con tribute to 10-20 % of the total energy output in the far- 
infrared (Pop e et al.l200 8). In Table[T]and Table|2]we summarize 
the data about CO emission in these objects collected from the 
literature. We note that the line luminosity in these objects is ex- 
tremely high, even for the higher J transitions, but also given 
the high value of the star formation rates, the resulting does 
not change too much. In addition, these very luminous objects 
are also very rare and therefore do not contribute significantly to 
the correlation term, which is due mainly to the most abundant, 
fainter objects. 

5.2.1. Dependence on the spectral resolution 

As shown in Sec. [H the spectral resolution of the instrument 
is extremely relevant in determining the amplitude of this sig- 
nal. This is one of the most interesting results of this work: we 
demonstrate analytically (AppendixlAb that the amplitude of the 
correlation term is expected to grow for lower values of Av/vobs 
to some threshold value, beyond which it reaches a convergence 
level and stays constant for even lower Av/vobs-s. This can be 
easily seen in Fig. |6] for the CO (1-0) line: the amplitude of the 
correlation term gains more than one order of magnitude as the 
relative spectral width is decreased to lO""*, but shows practically 
no change if this value is further reduced to 10 On the other 
hand, the Poisson term grows linearly with decreasing values 
of Av/vobs; therefore, the relative importance of the two terms 
changes as a function of the spectral resolution. 

We should keep in mind, however, that the actual amplitude 
of the Poisson term strongly depends on the ability of the ob- 
serving instruments to isolate and remove the bright individual 
sources. The bulk of the power in the amplitude of the Poisson 
term, in fact, is generated by rare and bright sources. Such 
sources can be removed from the maps observing the same re- 
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Table 3. Contribution to the amplitude of the correlation and 
Poisson terms of the CO (1-0) line from different intervals of 
flux. The observing frequency is 30 GHz and the spectral reso- 
lution is IQ--'. 
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gion with high enough angular resolution. This has practically no 
effect on the amplitude of the correlation term, which is gener- 
ated by the abundant low-flux sources. This is clear from Table|3] 
where we compute the contribution to Cf and Cf of the CO (1- 
0) line at 30 GHz from different flux decades; with a 1-mJy flux 
cut-off, one can drastically decrease the Poisson amplitude of 
70%, while decreasing the correlation of only 10%. 

The amplitude of the correlation term as a function of the 
spectral resolution is explored in more detail in Fig. |7]for the 
same line. Here we plot the value of the correlation term at dif- 
ferent multipoles I (in the range I -10 - 8000, corresponding to 
angular scales 9 ^ 2.5° - 1'). The signal grows uniformly and 
then reaches a plateau around a value Av/vobs ~ 10-^. The po- 
sition of this convergence point depends both on the redshift of 
the line and on the multipole /. In general larger scales (smaller /) 
and higher redshift lines seem to converge earlier, i.e. for higher 
value of Av/vobs- 

This can be explained as follows. The C/'s for the correlation 
term can be expressed as an integral in the line-of-sight compo- 
nent of the Fourier mode (k^) over the halo power spectrum, 

X*"""" dk 
■^Ph(k^,k,)\W(k,)\\ (10) 

where W(k,) is the Fourier window function along the radial di- 
rection, related to the spectral resolution of the experiment. The 
symbol kj^ represents the Fourier model included in the plane of 
the sky, and it is related to the I angular multipole via kj^ = l/r, 
with r the comoving distance to the observed redshift. The win- 
dow function W{k^) is such that k, 0, W{kS) — > 1, but for 
k^ » 2;r/Lsi- (with L^i- = c//-' Av/vinsn- the length associated to 
the spectral resolution), then W(k^) — > 0. In practice, the effect of 
this window function is therefore to change the upper limit of the 
integral in Eq. ( fTOl i by A:,, = InlL^^. Let us assume that the halo 
power spectrum Ph{k) is proportional to the linear matter power 
spectrum Pm{k), with k = (k^ + k]_y^^. For our choice of / and 
observing frequency (or redshift), at = the power spectrum 
is already being evaluated at ^'s comparable to or larger than 
~ 0.01 /;-' Mpc for which it shows its maximum. That means 
that, as long as k^ « k±, Ph(k) ^ Ph{k±) and when k- becomes 
comparable to or larger than kj_, then Phik) 0, since the power 
spectrum drops very fast (oc k^) when probing scales that entered 
the horizon during the radiation-dominated era. On these scales 
there should be strong non-linear corrections to the power spec- 
trum, which do not, however, significantly change the scaling of 
the integral in Eq. ( fTOl i. This integral can be approximated by 
Ci oc ksT Ph(k±) Vob.s/(Av) Ph(k±) for large widths Av/vobs, un- 
til it reaches a plateau for lower values of Av/vobs- This plateau 
will be reached later (that is, for higher values of k^^ or lower 
values of the spectral width) for correspondingly higher values 
ki. (or I). 
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Fig. 9. The contribution of the radio sources to the power spec- 
trum at 33 GHz (red), compared with the CO (1-0) line (blue). 
Different flux cutoffs 5 cm are applied: 10 mJy (solid lines), 1 mJy 
(dashed lines) and 0.1 mJy (dotted lines). The signal due to ra- 
dio sources is very sensitive to this value, while the power in the 
CO is much less affected, being generated by sources with fluxes 
much lower than 5cut- 



Observing in narrow spectral channels increases the sensi- 
tivity and observing time requirements. In the future, however, 
it will be possible to combine CMB measurements carried out at 
similar frequencies but with different spectral resolutions. For in- 
stance, future Planck L FI's channel at 30 GHz will scan a similar 
spectral range to CBI dPadin et alj|2002l) or VSA dWatson et al.1 
2003), but with a broader frequency response: Av/vobs - 0.2 for 
Planck, but Av/vobs ^ 0.05 for VSA and Av/vobs ^ 0.03 for CBI. 
Therefore, by studying the difference in the power spectra of the 
common areas, it will already be possible to establish constraints 
on the CO abundance. At 30 GHz we would simultaneously ob- 
serve two different slices of the universe at z ~ 2.8 and z ~ 6.7. 
But at lower frequency we have the opportunity to observe a sin- 
gle slice (i.e. at z ~ 4.8 for 20 GHz). Unfortunately, these limits 
would only be imposed up to a maximum multipole of I ~ 500, 
due to the relatively poor angular resolution of LFI's 30 GHz 
channel. Future ground-based experiments with arrays of detec- 
tors should be much more sensitive and with high enough spec- 
tral resolution could make such observations possible. 

In Appendix IB] we discuss a simple technique for separat- 
ing the line signal from the continuum, using different spectral 
resolutions. 

5.2.2. Dependence on the observing channel 

In this section we explore the amplitude of the signal in differ- 
ent observing frequencies. We consider the range of frequency 
5 - 300 GHz, which will be covered by several future experi- 
ments like SKA, PLANCK, ACT, SPT, and ALMA. In Fig. |8] 
we plot the ratio of the correlation signal to the primordial CMB 
signal at the fixed multipole / = 2000: here it is very clear how 
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the low-frequency range (v $ 100 GHz) is dominated by the 
CO emission, while at a higher frequency, the signal from dust 
is much stronger and dominates the signal from the CII 158 fim 
line. The sum from the CO lines is computed by assuming that 
the cross-terms are negligible. This shows that in the frequency 
range 30 - 60 GHz, the contribution from the first two CO lines 
might be the most important extragalactic foreground for the ob- 
servations of the primordial angular fluctuations by present-day 
and future experiments. The region around the curve correspond- 
ing to the first transition encloses the range of possible values for 
the signal in this line. The upper limit of this curve is calibrated 
on the Antennae. This local merging has a very bright CO (1-0) 
and a star formation rate comparable to M82. This results in a 
higher value of the !R ratio and in a higher amplitude of the sig- 
nal. As a lower limit, we used the IRAS 14378 galaxy. Similar 
ranges of uncertainty can be drawn around the other lines, ac- 
cording to the limits presented in Fig.[T] 

The contribution from dust emission in star-forming galax- 
ies falls rapidly towards lower frequencies, where the CO 
signal seems to dominate. However, in this range, the fore- 
grounds from compact radio sources (including AGN, quasars, 
BL Lac objects, af terglows of gamma ray bursts, etc., see e.g. 
Ide Zotti et al.ll2005l) is very relevant. Unfortunately, the cluster- 
ing properties of such sources are very poorly known, though 
several authors showed that the Poisson power is likely to be 
largely dominant over the co rrelation one (Toffolatti et al. 1998; 
iGonzalez-Nuevo et al.ll2005i) . The Poisson power, however, can 
be strongly reduced if the bright sources are excised from the 
map according to some flux cutoff. We show this in Fig.|9l here 
w e consider the radio source counts at 33 GHz, from the model 
ofSi ZottietalJ 12005). Integrating the counts we can estimate 
the Poisson contribution to the power spectrum. We compute it 
for several values of the flux cutoff' and we find that the signal 
from the radio sources can be reduced to the level of the CO ( 1 - 
0) line emission if a 5 cut - 0.1 - 1 mJy is applied. The same 
cut does not aflFect the CO contribution significantly, since at this 
frequency both its Poisson and correlation term are generated by 
much weaker sources. There are two important caveats that have 
to be considered here. First, the main limitation to such clean- 
ing procedure is the confusion noise. To reduce the radio sig- 
nal to the CO level, we need to apply a flux cuto ff" on the order 
of 0. 1 - 1 mJy. According to the counts curve of Ide Zotti et al.l 
(I2OO 5I). this corresponds to a density of ~ 20 - 200 sources 
per square degree. Therefore the confusion noise should not be 
an issue, even at very low ffuxes. Second, as mentioned above, 
the correlation term of the radio sources is considered negligi- 
ble compared to the Poisson power This is true, however, if the 
sources are n ot subtracted down to very low values of the ffux 
jToffolatti et a l. 1998). In this case the contribution of correlation 
from weak, unresolved objects might become important, and it 
is necessary to have an accurate modeling of their clustering. 

The galactic foregrounds are also independent of the spec- 
tral resolution, so they can in principle be separated from the 
CO emission fines. Moreover, their angular power spectrum de- 
creases rapidly towa rds the higher multipoles discussed here 
jTegmark et al.ll2000l) . especially at high galactic latitudes. The 
contribution from rotating dust is also characterized by a con- 
tinuum spectrum (iDraine & Laza rian 199 8i) and shows a peak 
around 20-30 GHz in our galaxy (.de Oliveira-Costa et al.l20()2l: 
iFinkbeiner et al.ll2004l) . which will be shifted towards lower fre- 
quencies in higher-redshift objects. 

To conclude this section, we recompute the correlation term 
for the CO (1-0) line using an alternative approach, based on 
a simple Press-Schechter distribution of the star-forming halos 



emitting in CO lines. This allows us to obtain a parallel con- 
strain to compare with the estimates of our merging model, in 
particular at the lower frequencies corresponding to the higher 
redshifts. 

Let us assume first that the star formation rate in halos scale with 
their mass 



= C(z) ■ M, 



(11) 



and the factor C{z) is obtained by normaUzing to the cosmic star 
formation rate density p*{z) 



r dn 

dM 

J dM 



C(z)M. 



(12) 



The integration must be limited to the halos that are actu- 
ally able to form st a rs. W e use the fit of proposed by 
iHernquist & Springell (l2003h . which uses a set of cosmological 
simulations up to very high redshifts. These simulation were per- 
formed with a value of cr°''* - 0.9, which is noticeably higher 
than the value adopted in this paper (cr^'^* - 0.74). Therefore, 
before computing the normalization in Eq. (fT2l i we rescale their 
fit taking this difference into account. Starting from Eq. (19) in 
Hernquist & Springel (2003), this rescaling can be written as 



p(trf,z) cr° 



p(o-"™,z) o-"' 



■expi-i 



D{z)cr\ 



.old 



1 - 



old \ 



cr 



(13) 



where 0-4 = criMA) and M4 is the mass corresponding to a virial 
temperature of lO'* K. This fit is accurate to better than 20 - 
25% for the redshifts corresponding to the peak of the signal 
and becomes even more precise at higher redshift. In either case, 
the error is much less than the other uncertainties of the model. 
Fig. E] shows the result for the first CO transition obtained using 
this simple model: the estimate is a factor of a few below the 
one obtained with the merging approach. The dependence of this 
normalization on the actual value of erg might open interesting 
possibilities of estimating the value of such a parameter from the 
amplitude of the signal. 

6. Angular fluctuations from fine-structure lines of 
atom and ions 

As mentioned in the introduction, the first generation of stars 
enriched the interstellar medium with significant amounts of 
metals, mainly carbon, oxygen, and nitrogen. The fine-structure 
transitions of these species are characterized by higher resonant 
frequencies with respect to the CO molecule. Therefore, their 
emission lines are more relevant for the high-frequency instru- 
ments, like PLANCK H FI, SPT, ACT, and ALMA. 

In an earlier paper, ISuginohara et al.l (Il999h presented es- 
timations of the collisional emission in fine structure lines of 
these species. Assuming a fraction of 10"^ the solar abundance 
of those species, and assuming that roughly half of them are 
in environments with super-critical densities, they find that this 
emission should introduce an average distortion at the level of 
10"^, whereas on scales of 1 /i^'Mpc, the relative intensity fluc- 
tuations should be a factor of 10 larger. 

In this section we compute the signal due to the clustering 
of the sources responsible for the line emission. We repeat the 
approach used for the CO lines. To calibrate the "R ratio for the 
fine-structure transitions of atoms and ions, we use the results 
of Malhotra et al. (2001), who observed 60 star-forming galax- 
ies with ISO-LWS. They especially selected the objects with a 
smaller angular size than the instrument beam, therefore the line 
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Fig. 10. The correlation signal of the most significant atomic 
emission lines, for a spectral resolution Av/vobs - 0.33 at 
353 GHz, corresponding to one of the PLANCK HFl's detectors. 
The black line is the primordial signal of the CMB, the orange 
line the signal from dusty merging star-forming galaxies. Other 
colors identify the lines as indicated by the labels. See Table |4] 
for details. 
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Fig. 11. The ratio of the correlation signal from atomic and ionic 
line emission to the primordial CMB signal at I = 2000 as a 
function of frequency and for a spectral resolution Av/vobs - 
10"^. The orange line is the signal from dust emission in merging 
galaxies. The Cll line clearly dominates the other species, but is 
below the dust signal even at this spectral resolution. 



Table 4. Average 'R for the 7 transitions of the sample from 
iMalhotraetanOOOlh . 



Species 


Wavelength [/im] 


Average K 


CII 


158 


6.0 X 10* 


01 


145 


3.3 X 10^ 


Nil 


122 


7.9 X 10^ 


OIII 


88 


2.3 X 10* 


OI 


63 


3.8 X 10* 


NIII 


57 


2.4 X 10* 


OIII 


52 


3.0 X 10* 



fluxes presented in their paper refers to the total line emission 
from such galaxies. As for the low-redshift sample of CO lines, 
we obtain the FIR luminosities of these objects fr om the IRAS 
revised bright galaxies sample { S anders et alJ2003h . For each of 
the 7 species presented in their work, we computed the geomet- 
ric average of the R ratio in the sample as 



R^ 



l/N 



(14) 



The resulting numbers are listed in Table |4] 

In Fig. [To] we show the correlation signal of these lines at 
353 GHz and for a spectral resolution Av/vobs = 0.33, close 
to the performance of one of PLANCK HFl's channels. In 
the same figure we also show the correlation term from star- 
forming merging d usty galaxies computed in our previous paper 
jRighi et al.ll2008l) : it is clear that this signal is still the domi- 
nant one. Also in this case, however, it is possible to increase 
the amplitude of the correlation term by using very good spec- 
tral resolutions. We show this in Fig. (TT] for a spectral resolu- 
tion Av/Vobs - 10"^: the relative contribution of the Cll 158 //m 
line increases with frequency when compared with the contin- 
uum emission from dust in merging star-forming galaxies. Other 
lines, plotted in the same figure, have a much lower amplitude 
and are therefore neghgible. 



6.1. Comparison with similar works 

An estimate of th e fluctuations arisi i ig from atomic line emission 
was presented bv lSuginohara et aT I (119991) . When looking at the 
Cll line at z = 10, we find that our results predict emission am- 
plitudes that are significantly (a factor of ~ 10^) less than theirs. 
Most of this disagreement is due to the higher value adopted for 
the density of ions emitting in the line (which depends on the 
assumed metallicity and ionization/excitation fractions), while 
having a different set of cosmological parameters introduces an 
offset of a factor of a few. In our work, the actual abundance of 
ions in the upper state is provided by the calibration of our model 
wit h observational da ta. 

iBasu et all (12004 ') analyzed the effect on the CMB primordial 
spectrum induced by the resonant scattering of heavy elements in 
the intergalactic gas, expelled by supernova explosions, galactic 
winds, and jets. They demonstrate that this effect is the sum of a 
damping of the original fluctuations plus the generation of new 
anisotropy. They obtained a very simple analytical expression 
for the change induced in the angular power spectrum on small 
scales 



6Ci 



"2Ts(;atC;; 



(15) 



where Tscat is the scattering optical depth of the line. The linear 
scaling of dCi with the optical depth produces quite strong effect 
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Fig. 12. The effect of the scattering of the CO (3-2) line at 
70 GHz (squares), according to the model of iBasu et al. 1 (12004 . 
The solid line is the correlation term from the emission in the 
same transition, for a spectral resolution Av/vobs = 10"^. 

even for a low value of r. In Fig.[T2]we plot the 6Ci for the third 
transition of CO and compare it with the correlation signal of the 
emission from the same line. For 1 % solar abundance, the scat- 
tering optical depth of this line at 70 GHz is T3_2 = 1 .94x 10 It 
therefore introduces a change in the primordial power spectrum 
on the order of ~ 4 x 10""*. This term, however, dominates on 
large scales, while it drops towards higher multipoles, following 
the decrease in the primordial spectrum, where the correlation 
signal dominates. 

More recently, iHernandez-Monteagudo et alj (l2007l |2008|) 
considered the pumping effect of the ultraviolet background on 
the OI 63 jim line and the resulting distortion on the CMB. They 
compute the correlation signal associated with the clustering of 
the first star-forming objects. The amplitude they find is on the 
level of 10"* (juK)^, and it is therefore lower than the signal de- 
scribed here. However, this effect would show up at higher fre- 
quencies (v ~ 400 - 700 GHz), well within the spectral coverage 
of ALMA, whose sensitivity should be able to put constraints on 
this effect. We notice that the physical environments where both 
the coUisional emission and the UV-induced emission take place 
are very similar, and therefore both effects constitute different 
probes for the same scenarios. 

7. Conclusions 

The same star-forming activity that causes reionization of the 
IGM at z ~ 10 leaves an imprint on the CMB by means of several 
mechanisms (e.g., Thomson scattering on the ionized gas, IR 
emission from dust particles that reprocess UV radiation, reso- 
nant scattering on metals produced by the first stars, etc.). In this 
work, we have focused on the impact that emission on molecu- 
lar and atomic and ionic lines have on the angular power spec- 
trum of the CMB. We put particular emphasis on the emission 
of CO rotational lines, since their effect is particularly strong in 



the 20 - 60 GHz frequency range, although we have also consid- 
ered the collisional emission of lines corresponding to species 
Uke CII, OI, Nil, and OIII. 

If the star formation activity follows the halo merging his- 
tory, then the spatial distribution of molecules and metals should 
closely resemble that of the halos where they were produced; 
therefore the emission pattern generated by these molecules and 
metals should show a similar clustering pattern. The advantage 
of looking at the emission of particular lines is that each observ- 
ing frequency probes a given redshift shell, whose width is deter- 
mined by the experimental spectral resolution (Az/z ~ Av/vobs). 
The anisotropy introduced by the clustering of the sources will 
be optimally measured if both the angular and the spectral reso- 
lutions of the observing experiment are able to spatially resolve 
those scales corresponding to the clustering (~ 15 - 25 /i"' Mpc, 
in comoving units). On top of these clustering-induced fluctua- 
tions, we have to add the fluctuations associated to the Poisson 
nature of the sources, which contribute down to much smaller 
scales (the typical source size) and for which further improve- 
ments on Av/v result in a larger amount of measured anisotropy. 
In this context, all other signals that contribute in similar spectral 
ranges (dust IR emission, synchrotron radiation, intrinsic CMB, 
etc.) show a very different behavior with changing spectral reso- 
lution, because they merely remain constant. This will provide a 
powerful tool for distinguishing between the fluctuations charac- 
terized in this paper from all the rest and an efficient way to per- 
form tomography of reionization at different frequency bands. 

Due to the large number of uncertainties when building 
a theoretical model for the collisional emission on molecular, 
atomic, and ionic lines, we used existing measurements at low 
and high redshift to calibrate our model. Our approach is simi- 
lar to lRighi et al.l (^008), where our halo merging history model 
assigns a value of the the star formation rate to a given object 
and where this rate is assumed to follow a linear scaling with 
the luminosity under study (in this case, the luminosity in the 
various lines considered). For the CO rotational lines, we col- 
lected a sample of local and high-redshift objects of different 
sorts (including local mergings, bright infrared sources, AGN, 
QSOs, and radiogalaxies) to calibrate the 'R = L/M ratio. We 
have considered, for each transition, a broad range of possible 
values to give an upper and lower limit estimate of the expected 
amplitude of the fluctuations. 

Our results show that CO emission is already stronger than 
dust emission in the frequency range 10-70 GHz if Ay/v = 10"^. 
Moreover, if one makes use of the different behavior of the 
CO signal with the observing spectral resolution, then it should 
be possible to project out continuum signals like CMB, radio, 
and dust emission with greater confidence. In particular, the fre- 
quency interval ranging from 30 up to 50 GHz is of special in- 
terest for the first three CO rotational lines, since they probe the 
redshift range z ^ 6 - 10, crucial for our understanding of reion- 
ization. Atomic and ionic emission lines are more important at 
higher frequency, but only the CII 158 /urn gives a non-negligible 
contribution to the amplitude. 

This brings us to the conclusion that the emission on CO 
and CII lines studied in this paper provides a new window into 
reionization, constituting a bridge between the low-frequency 
observations pursuing the HI 21 cm fluctuations in the radio 
range, and the high-frequency observations targeting fine struc- 
ture l i nes in metals and ions like CII, OI , or OI II jBasu et alJ 
12004 IHernandez-Mont eagudo eTal] l2007l 12008'). These three 
windows should allow tomography of the same cosmological 
epoch to be performed, but should be affected, in general, by 
very different contaminants and systematics. Therefore, the com- 
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bination of these three perspectives should provide a consistent 
picture of an epoch that, to date, has remained hidden to obser- 
vations. 
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Appendix A: The induced angular power spectrum 

A. 1 . The measured intensity 

In this section we compute the angular power spectrum induced 
on the CMB by the emission of these lines. We follow the 
line of si ght approach (los) already outlined in Appendix A of 
iRighi et al. (2008) and in Hernandez-Monteagudo et al. (2008!), 
although here we place more emphasis on particular aspects re- 
lated to the observing frequency resolution. Let B(h, h') be the 
PSF of the observing beam on the sphere unit vector n' while 
pointing along n. Let Vobs also be the observing frequency of the 
experiment that measures the emission of a line of rest frequency 
V and therefore probes the universe at a redshift 1 -i- z = v/vobs- 
The spectral intensity introduced by such line reads 

Aly - J dr J dn B(h,n') J dv' (f>yi„f,tr x 



dLv', 



(A.1) 



v'(l+z) 
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Here, 0v',insti is the frequency response of the experiment, 
and dn/dLy>(i+.) the luminosity function of halos hosting the 
metal/ionic species responsible for the line emission under study. 
The emissivity on the line at the center of the halo is given by 
7v'(i+z')' which is diluted by cubic power of the scale factor a. The 
function W{y - r') provides the profile of the density distribution 
of the species in the halo. The vector r' is determined by n' and 
the frequency v' providing the effective distance to the observer: 
r' = r'(v') n' . At this point, al l spatial coordinates are physical. 
As shown inl Righi et alj (1200 8^. the convolution in Eq. dA. ll ) can 
be rewritten in terms of the source luminosity 



My = J dr JdnB(h,A') J dv' (py', 
dh(r') iv'(i+z) 



dLyi(i+r) 



dLy' 



An 



(A.2) 



The new number density dhldLyni+r) is just the previous one 
smoothed on scales corresponding to the typical halo size. The 
spectral luminosity can be expressed as a product of the bolo- 
metric luminosity on the line and the emission spectral profile, 
Ly'iUz) - hxA^v'ix+z)- We approximate this profile as a top hat 
function of the thermal width of the line, (i/'v'(i+,) = l/(Av)th 
if |y'(l + z) - v| < (Av)th, lAv'Ci+z) - otherwise), and assume 
that this width will be much smaller than the instrumental one 
((Av)th/v «: (Av)instr/vobs)- This enables the replacement of the 
integral along r by c//"'(z)(Av),h/v, yielding 



My = ^ dnB{n,h') J dv' 



'v^instr 



dL 



v'(\+z) 



dn(r') iv'(i+z) 



dLyl 



V(l+z) 



4-7T 



cH-Hz), 



where Ly = L\,o\/v. This cancels any dependence on (Av)th. Note 
that the halo number density is evaluated at r' and that the in- 
tegral along the line of sight is carried out under the frequency 
response of the instrument: 

My = J dr"P{r-r') J dh B(n,h') x 
dn(r') iv'(i+z: 



dL 



'y'(l+z) 



dLyi(i+,) An 



' cH-'{z). 



(A.4) 



The product fir - r') = dv' /{dr') 0v',insti- constitutes the radial 
profile of the los integration. We next attempt to show the im- 
pact that the shape of the radial (and angular) responses of the 
experiment have on the measured intensity. If we convert n' into 
the transversal spatial component by introducing r'^, Eq. jAAj 
reads like a convolution 



A/,(r) = j dr' Sir - r') /v'(i+z)(r'), 
with S(r - r') = !P(r - r')/r'^ and 

£v'(i+z) dh{r') 



(A.5) 



/,,(i+,)(r') = c//-'(z) 



J dLy'(l+,) 



An dL 



(A.6) 



y'(l+z) 



Here, S(r - r') represents the 3D window function of the in- 
strumental response and is responsible for some smoothing of 
the original signal 7y'(i+2)(r'). This becomes more obvious if we 
rewrite the convolution in Fourier space 



A/v(r) 



r dk 
"J (2^ 



- exp(-/k-r)%/v-(i+j)(k). 



+z) 



(A.8) 



If the instrument is not sensitive to small scales (due to a poor 
angular or spectral resolution), then at large k the will van- 
ish, suppressing the contribution of 7v'(i+,)(k) at the same k's. 
This is also reflected in the difference scaling of the power spec- 
tra versus the spectral resolution of the experiment. Indeed, if 
we compute the power spectrum of A/y(r) in terms of the halo 
power spectrum, we obtain 

• \ 47r dLy'i^i 

with dh/dLy the average number density of halos per unit lu- 
minosity. If the halos are Poisson-distributed, then Phik) oc 
l/{dh/dLy>^]+,y) up to very large k-s, and the high-fc integral is 
actually limited by the window function Sy^ (as long as k corre- 
sponds to scales larger than the source size). In contrast, if halos 
are clustered so that their power spectrum is proportional to the 
linear matter density power spectrum Pm(k), then we find that 
Pm{k) — » for some large k, and the increase in spectral/angular 
resolution of the experiment will not make any difference at this 
high k range. Therefore, improving the spectral resolution helps 
for increasing the Poisson contribution, but does not change 
the contribution from the clustering term beyond some large 
kc — 2n/Lc, with Lc the typical clustering scale. 



A.2. The angular power spectrum 



Our starting point is Eq. iA.Ai above, which can easily be con- 
verted in comoving units by absorbing the a^ factor within the 
definition of comoving number density. An a factor must then 
be added for the line element dr, but this cancels out with the 
(A. 3) ' foi" the profile 'P. To compute the angular power spectrum, 
we have to characterize the clustering properties of the emitting 
halos. We recall that 



dn r dn ^ 

^^-jdM-GiM,Ly,,^,), 



(A.9) 



where the function G(M, Lv'(i+,)) provides the fraction of halos 
present in the mass interval [M, M + dM] (given by the mass 
function dn/dM) that have recently experienced a major merger 
(see Righi et al. ( 2008)), and hence given rise to significant emis- 
sion in the line(s) of interest. The spatial correlation function of 
halos in the mass range [M, M + dM], n^iM, x) reads 

dn dn 
{nh{MuXi)nk{M2,X2)) = — (Mi, jci) — (M2, ^2) 
dM dM 



dn , 
+ —(MuXx)S',j{xi 



■ xi) 6d(Mi - M2) 



dn dn 
+ ■777(^i'-^i):7T7(^2,X2) 
dM dM 

X biMuZ[xi])b(M2,z[X2])U(Xi -X2). 



(A. 10) 



The symbol 6d stands for Dirac delta and b(M, z) is the mass and 
redshift dependent bias factor that rela tes the halo and the matter 
linear correlation function dMo & W hite 2002). The first term 
on the righthand side accounts for the Poissonian fluctuations in 
the number counts, whereas the second term accounts for the 
dependence of the halo number density on the envir onment- 
Ha ving this present, just as in Appendix A of Rig hi et aP 
(I2OO8I) . it is possible to write the angular correlation function of 
the intensity fluctuations as 
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with P/(ni ■ Aa) the Legendre polynomia of order I and where 
and cf are the Poissonian and the correlation terms of the 
/-th angular power spectrum multipole. Then, Cj" can easily be 
found to be 



(Av)„ 



Vobs 



liLv'd+z) cH (z) 



'v'(l+z) 



dn 



4-n I dLy'(i+z) 



(A. 12) 



On the other hand, the correlation term can be expressed in 
terms of a A;-space integral of the initial scalar metric power 
spectrum P,:,(k) times a sq uared transfer function, just as in 
ISehak & Zaldarriagal d 19961) . 



C 



f = ^ J k^dkP^(k)\A,(k)\\ 



with the transfer function A/(A:) given by 
Ai(k)^ j drjiikr) P(r) [S(r)6k]. 



(A. 13) 



(A. 14) 



In this equation, ji(x) is the spherical Bessel function of or- 
der /, P{r) is the instrumental profile function as defined after 
Eq. iAAi . 6k is the fc-mode of the dark matter density contrast, 
and the function S (r) is defined by 

Sir) s J dLv'(i+.)c/MG(M,Lv'(i+,)) 

^ 777 — r^c//-'(z)MM,z[r]). (A.15) 
dM An 



An optimal estimate for the line-emission amplitude can be 
obtained by minimizing the quantity 



ijjm 

[di^-&(m)-m,))], 



(B.3) 



i.e., by finding the value of a that satisfies dx^ Ida - 0. The co- 
variance matrix of the combined channel noise is given by Nijjm- 
The minimization yields 



" ~ lluM {m) - fiSj)) Nr.^, {f(6i) - m„)) ' 
with a formal error 

1 



" Y.ij,i,n {m) - fisj)) NT^,^ im) - m„y) ' 



(B.4) 



(B.5) 



Appendix B: Observational method 

In this Appendix we propose a simple method for extracting the 
line-induced signal from a continuum at a given frequency, by 
the use of observations at different spectral resolutions. 

Our model for the measured signal (in either real or Fourier 
space) is given by an array of channels at different spectral 
resolutions (5/ = (Ay),7v, / = 1,A^): 

~Si^af{6i) + C + Ni. (B.l) 

The vector Si should sweep the relevant spectral resolution range 
shown in Fig.|7] (5, € [0.1,5 x 10"^]). The value of a provides 
the amplitude of the line-induced emission, whereas the function 
f{5i) contains its dependence with respect to the spectral resolu- 
tion. A priori, this function can be characterized if the redshift of 
the line-induced emission is known. Signal C refers to all com- 
ponents that do not show an explicit dependence of spectral res- 
olution (e.g., intrinsic CMB, dust, synchrotron, etc.) and whose 
value should remain constant (at least down to a few percent) 
in every channel. Finally, A^, refers to the instrumental noise in 
each channel, and its covariance matrix is assumed to be char- 
acterized. With this in hand, it is convenient to take differences 
between measurements in different channels having distinct val- 
ues for /(v), 

dij = ~Si -Sj^a {f(6i) - f(6j)) + Nij, (B.2) 

where Nij denotes a combination of the noise from both chan- 
nels. 



